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Abstract 

The  density  of  states  in  delta-profiled  2D  and  ID  quantum  systems  is 
calculated.  It  is  shown  that  there  are  smooth  crossovers  in  the  density  of 
states  from  a  3D  square-root  behavior  to  a  2D  steplike  behavior,  and  from  a  2D 
to  a  ID  sawtooth-llke  behavior,  as  the  confinements  Increase. 


PACS  Numbers: 


73.40.Gk,  71.20. -b 
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I.  Introduction 

Quantum- we  11  structures  are  being  extensively  studied  for  applications 
in  ultra-large-scale  integrated  circuits  and  high-speed  optoelectronics.  As 
for  these  systems,  the  lateral  confinement  of  originally  quasi- two-dimensional 
(Q2D)  electron  layers  to  submicron  dimensions  has  made  possible  the 
realization  of  quasi -one-dimensional  (QlD)  electron  systems[l].  Among  them, 
resonant- tunneling  structures  are  of  great  interest  not  only  because  of  their 
potential  applications,  but  also  for  the  underlying  basic  physics.  A  number 
of  theoretical  studies  have  been  carried  out  recently[2-5] .  For  example,  one 
of  the  authors  has  calculated  the  density  of  states  and  dwell  times(3,4],  and 
Bahder  et  al  have  calculated  the  local  density  of  states  for  a  simplified 
model  with  quite  rigorous  results [5], 

Despite  extensive  investigations  of  quantum-well  structures  to  date,  we 
are  not  aware  of  any  study  of  the  cross  sectional  local  density  of  states  of  a 
quantvuB-well  wire  structure.  To  solve  this  problem,  we  consider  an  artificial 
structure,  the  so-called  two-direction  double-barrier  resonant-tunneling 
structure.  In  Section  II  we  evaluate  the  eigenfunctions  and  eigenvalues  from 
an  effective-mass  Schrodinger  equation.  For  the  delta-profiled  potential 
model  using  these  results  in  Section  III,  the  local  density  of  states  has  been 
determined.  In  Section  IV,  integrating  over  the  well  volume,  we  calculate  the 
density  of  states  in  the  well,  where  we  show  crossovers  of  the  density  of 
states  from  3D  to  ID  via  2D. 


II .  Theoretical  model 

Ve  consider  a  typical  two-direction  double -barrier  resonant- tunneling 
structure  (DBRTS)  consisting  of  two  thin  (-  50  A)  Al^Ga^^  ^As  layers,  separated 
by  a  thin  GaAs  layer  in  both  directions.  The  potential  is  expressed  by 


Uly  t 

ir' 


. y  CciOS 
R'K’/or 
ipeolel 


□  □ 
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V(y,x)  -  V  (y)  +  V  (z) 

y  2 

-  V^{5(y+b)  +  5(y-b))  +  V^{Siz+a)  +  S(z-a))  (1) 


In  this  expression,  the  four  Al^Ga^  ^As  potential  barriers  have  been  replaced 
by  S-functions  with  strengths  and  in  the  y-  and  z-direction, 
respectively.  The  parameter  (i  -  0  or  1)  is  given  by 


V.  -  d.AV  . 

1  1  Cl 


(2) 


where  d.  are  the  barrier  widths  and  AV  .  are  the  c o nduc  t  i  o  n  -  b and 

1  Cl 

discontinuities.  This  additive  potential  form  corresponds  to  a  rectangular 
quantvun  wire  with  cross  section  a  x  b  when  and  are  very  large  and  a  and 
b  are  less  than  the  deBroglie  wavelength  (A^)  of  the  electron.  The  quantum 
wire  with  circular  cross  section  also  has  a  quasi-one-dimens ional  character, 
but  the  wavefunction  in  the  confinement  direction  is  the  Bessel  function. 

We  now  solve  the  time- independend  Schrddomger  equation  with  the 
Hamiltonian 


m2  0 

H  -  -  5^  +  V(y.z) 


(3) 


where  m  is  the  effective  electron  mass  at  the  bottom  of  the  GaAs  conduction 
c 

band.  In  order  to  deal  with  a  finite  density  of  states,  we  must  take  our 
structure  within  a  large.  Impenetrable  rigid  box  extended,  say,  from  -L/2  to 
L/2.  With  these  boundary  conditions,  the  Schrddinger  equation  is  separable 
with  the  additive  potential  form,  and  then  we  can  write  the  wavefunction  in 
the  product  form 
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♦(r)  -  L 


exp(ik^x)  'Ky)'l'(z) 


(4) 


where  -  2ffn^/L  and  takes  the  integer  values  0,  1,  2,  ...  .  The  y-  and 
z-parts  of  the  wavefunction,  ’l'(y)  and  ^'(z),  satisfy  the  reduced  equations 


2m 

*"(7)  +  ^  [E^-V(y)],6(z)  -  0  (5) 


2m 

♦"(z)  +  -f  [E  -V(z)]^(z)  -  0  ,  (6) 


where 


E  +  E  -  E  -  )f(\^/2m 
y  z  X  c 


(7) 


Here,  E  Is  the  total  energy  corresponding  to  the  Hamiltonian  H,  and 
E  (E  )  is  the  energy  eigenvalue  of  Eq.  (5)  ((6)).  Because  of  the  symmetry  of 

y  ^ 

our  system,  it  is  convenient  to  write  the  wavefunction  in  terms  of  even  and 
odd  functions  as 


(y) 


Aj^(ky)cos(kyy)  , 

A2(ky)cos(k^)  +  Aj(ky)sin(k^)  , 


0  <  y  <  b 

(8) 


b  <  y  <  L/2 


and 
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B^(k^)sin(k^) 


0  <  y  <  b 


^k  - 


(9) 


B2(k^)cos(k^)  +  B2(ky)sin(k^) 


b  <  y  <  L/2  . 


The  z- components  of  the  wavefunctions  can  be  written  in  similar  forms  as 


C^(k^)cos(k^z)  , 


0  <  z  <  a 


♦ek  - 

z 


(10) 


C-(k  )cos(k  z)  +  C_(k  )sin(k  z)  ,  a  <  z  <  L/2 

2  z  z  i  z  z 


and 


0  <  z  <  a 

(11) 

a  <  z  <  L/2 

Uhen  we  apply  the  boundary  conditions  to  the  y-components  of  the  wavefunction, 
we  get  the  equations  for  the  bound  states  of  even  and  odd  parity, 
respectively,  as 


♦ok  “ 

z 


D-(k  )sin(k  z)  , 
i  z  z 


D2(kz)cos(kzZ)  +  D2(k^)sin(k^z)  , 


h 


cos(kyb)sin(kyL/2-k^b)  +  cos(k^L/2)  - 


(12) 


and 
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^  sin(k  b)sin(k  L/2-k  b)  +  sin(k  L/2)  -  0 

y  y  y  y 


(13) 


where  7,  -  2m  V, /W  .  The  coefficients  A.  and  B.  are 
1  c  1  1  1 


^3  2 

7^  -  ^  cos^(k  b) 
A,  k  V 

1  y 


(14) 


1  -  ^  sin(2k  b) 

y 


(15) 


1  +  —  sin(2k  b) 

y 


(16) 


^2  ^1  2 
o  sin^(k  b) 

B-  k  y 

I  y  ^ 


(17) 


With  the  same  calculation  for  Eq.  (6),  we  find  similar  results  for  the 

conditions  of  the  bound  states  and  coefficients  and  by  replacing  b  with 

2 

a,  ky  with  k^,  and  7^  with  7q.  which  is  2m^V^/>(  .  From  the  normalization 
conditions,  A^^,  B^,  and  can  be  determined  as  shown  for  A^  and  B^^  in  Ref. 
[5],  The  energy  eigenvalues  corresponding  to  Eqs.  (5)  and  (6)  are  given  by 


(18) 


where  a  (*-  e  or  o)  labels  the  state's  parity. 

Taking  7^^  and  7^  both  to  be  equal  to  zero,  which  is  identical  to  the 
limit  where  the  j-functions  are  placed  on  the  boundaries,  that  is,  b  -  a  - 
L/2,  one  recovers  A^/Aj^  -  -  ^3^^!  “  ®2^®1  ”  ^2^^1  ”  ”  ^2^^^!  ” 
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Dj/D^  -  1  and  -  8^  -  -  [2/L]  which  is  the  result  for  the  motion 

of  a  3D  particle  in  a  rigid  box. 


III.  Local  density  of  states 

The  local  density  of  states  (DOS)  in  the  DBRTS  has  been  obtained  in 
various  cases  [6] .  It  can  be  defined  in  the  two-direction  case  as 


N(y,z;E)  -  -  ~rniG(r,r:E) 

If 


k  a/3  k  k 
X  y  z 


Iv I\k  ■ 

y  z 


where  the  factor  of  2  implies  spin  degeneracy,  G(r,r':E)  is  the  single¬ 
particle  Green's  function,  and  a  and  /3  (-  e  or  o)  label  a  state's  partiy. 
When  the  system  size  goes  to  infinity,  we  can  change  the  summation  into  the 
appropriate  integration  because  the  density  of  allowed  wavevectors  becomes 
2jr/L: 


with 


5  “  2^^1  ®1  (A2.Bj)cos(2kyy)]  , 


}  “  2^^!  °1  (Ci-D^)cos(2k^z)  . 
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Integrating  Eq.  (20)  over  k^,  ve  get 


N(y,z;E) 


[A^  (A^-B^)cos(2k  y) 


X 


[cj  +  D?  +  (cJ-Dj)cos(2k  z)] 
—1 -i - 1 — 1 - S 

[E  -  [K^Sni  ](kW)]'’ 

t  y  ^ 


(23) 


The  wavefunction  coefficients  A^,  and  are  given  by 


lim  (L/2)A^p/b) 
2b/L-*0 


2 

~2 - 2 - 2 - ^ -  ■  ^e^P^ 

p  +  U^cos  (p)  -  p  sin(2p) 


(24a) 


11m  (L/2)B^(p/b)  -  — - - - ; - 2 -  .  p  ,p)  (24b) 

2b/L-»0  ^  +  U^sin(p)  p  sin(2p) 


lim  (L/2)cJ(q/a) 
2b/L-*0  ^ 


2 

_ a _ 

q^  +  U^cos^(q)  -  q  sin(2q) 


Gg(q) 


(24c) 


11m  (L/2)D?(q/a) 
2a/L-*0  ^ 


2 

- a _ 

2  2  ^ 

q  +  U  sin  (q)  -  U 
^  o  ^  o 


q  sin(2q) 


Go(q) 


(24d) 


where  p  -  k  b,  q  -  k  a,  U,  -  y.b  and  U  -  7  a.  Furthermore,  when  we  allow  V- 
y  ’z  l  l  0  0  1 

or  Vq  to  go  to  zero,  N(y,z';E)  becomes  N(z;E)  or  N(y;E),  which  is  that  of  the 
one-direction  DBRTS. 

Equation  (23)  shows  the  cross-sectional  local  DOS  of  the  quanttim-well 
wire,  the  so-called  two-direction  DBRTS,  which  consists  of  two  parts,  each 
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coming  from  even  and  odd  parity,  respectively.  In  both  the  y-  and  z- 
directions,  it  shows  sinusoidal  behavior.  In  the  limit  of  -►  0  and  Vq  -*  0, 
the  functions  F(p)  and  G(q)  are  unity  and  N(y,z;E)  becomes  the  DOS  of  a  free 
electron  in  a  box  of  volume  uabL. 


IV.  Crossovers  of  the  density  of  states 

Ue  now  consider  the  DOS  in  the  well,  N(E),  which  can  be  calculated  by 
taking  the  integral  over  the  well  volume. 

|•L/2  pb  pa 

N(E)  -  8  dx  dv  dz  N(y,z;E)  (25) 

The  result  is 


L2a  pbrr 


N(E)  - 


Jo 


pbfj 

r  -p  If. 


(P)  +  F^(p)  +  (F^(p)-F^(p))stn(2p)/2p] 


pa|f>^-(p/b)^)^  G^(q)  +  G^(q)  +  (G^(q) -G^(q)  )sin(2q)/2q, 

dq  [-S - 2- - ^ - J  ^  (26) 

'O  {i,^  -  (q/a)^  -(pA)  1 


where  -  2m^E/)(^. 

Lac  us  evaluate  the  DOS  for  a  few  extreme  cases. 

(1)  30  case 

This  corresponds  to  the  limits  «  1  and  «  1,  i.e.,  F^(p)  ”  F^(P) 
-  G^(q)  ”  ^***  arrive  at 


8m  L  pbfj  fa(ij^-(pA)^)^  j  ,  5  u 

N(E)  -  -fr  dp  dq  I/I7  -(q/a)^-(pA)  ) 

w  K  -*0  -^0 
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.  2m 

4abL  — f  e’^ 
In  iT 


which  is  the  well-known  DOS  of  a  3D  free-electron  gas  with  volume  4abL. 

(ii)  2D  case 

In  this  case,  either  or  goes  to  zero  while  the  other  goes  to 

infinity,  i.e.  ,  U  «  and  U,  ->  0  or  vice  versa,  such  that  Eq.  (26)  becomes 
^  o  I 


*4Uil  pU 

N(E)  - 

K  ]ll  Jo 


4Lni^  pb{fr^-(q/a)^}‘’ 


dp  l/(n^-(p/b)^-(q/a)^)‘* 


X  dq  (G  (q)  +  G  (q)  +  (G  (q)  -  G  (q) )sin(2q)/2ql  (28) 

J  y\  6  O  ®  O 


For  U  -*  •,  Che  DOS  becomes 
o 


4Lb  b 


N(E)  - 


1 C  «l 


6{q- (n+l/2)»r  5{q- (n-i-l)jr)  +  [5{q- (n+l/2)jr) 


-  6|q-(n-Kl)ir)  1  [sin(2q)/2q] 


m  2bL 


IC 


dq{5{q-(n+l/2)s’)  +  S (q- (n+l))r)  ] 


m  2bL 


—  2  ’ 


2  2  2. 

where  9  is  the  unit  step  function  and  E  -  ir  h  /(8m  a  ).  Here,  we  have  used 

o  c 

the  fact  that  can  be  represented  in  terms  of  5 -functions  when  Uq  goes  to 

infinity  (see  Fig.  1). 

(iii)  ID  case 

This  corresponds  to  and  both  going  to  infinity  so  as  to  confine 
the  motion  of  the  electrons  to  just  one  direction: 


L2m  pbi;  ^ 

-  '"3^1  dp  JT  y  [5 {p- (m+l/2)ir)  +  (S{p- (m+l)iT)  ] 


?  fa{»,^-(p/b)^)l/2 


^  I  r 


'n+i/Z)ff I 


2  2  2  h 

[/  -  (q/ar  -  (9/b)  y 


r2m  t  h 


"  }  [l/[E.()<V2m^){((m+l/2)/b)^  +  ((n-H/2)/a)^)]* 

^  m,n 


+  l/[E-(xV/2m^){((m+l/2)/b)^  +  ((n+l)/a^}  ]** 


+  l/[E-(xV/2m^){((nri-l)/b)^  +  ((n+l/2)/a)^)  1^* 


+  l/[E-(xV/2m^)(((m+l)/b)2  +  ((n+l)/a) )  ]^)  ]**] 


We  have  again  used  the  fact  that  F^(p)  behaves  as  5-function  when  goes  to 


In  order  to  find  the  crossovers  of  the  DOS  in  the  well  from  3D  to  2D 


graphically,  we  modify  Eq.  (26)  as 
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2  2 
2bLm 


N(E) 


t/2(E/E^)  - 

dq  (G^(q)  +  G^(q)  +  ( G^(q) -G^(q) ) sin(2q)/2q ] 


(31) 


u 

where  jr/2(E/E^)  -  317  and  G^(q)  is  given  by  Eqs.  (24c)  and  (24d).  Similarly, 
the  equation  which  shows  that  the  crossovers  of  the  DOS  in  the  well  from  2D  to 
ID  is  expressed  as 


2Lm 

c 


«  t/2(E/E^)^ 

N(E)  -^1  dq[G^(q)  +  G^(q)  +  { G^(q) -G^ (q) ) sin( 2q)/2q] 

m— 0 


[(ir/2a)^E/E^  -  { (m+l/2)>rA)^  -  (q/a)^]‘* 


((ir/2a)^(E/E^)  -  {<m+l)x/b)^  -  (q/a)^]^*-* 


(32) 


Figure  2  shows  the  graphical  result  of  Eq.  (31),  namely,  the  crossovers 
of  the  DOS  in  the  well  from  3D  to  2D.  In  this  case  we  take  -  0,  changes 
from  0  to  20,  and  E/E^  varies  from  0  to  8.  On  the  other  hand.  Fig.  3  shows 
the  transition  of  the  DOS  in  the  well  from  2D  to  ID.  For  the  sake  of 
convenience,  we  take  a  -  b,  to  go  to  infinity,  and  to  vary  from  0  to  20. 
Higher  values  of  correspond  to  increased  sharp  peaks  of  the  DOS  of  the  ID, 
quantxim-wire  case.  In  this  case,  our  result  recovers  the  well-known  sawtooth 
type  DOS  diverging  at  values  of  E/E^  -  2,  5,  8,  10,  ...,  which  is  in  good 
agreement  with  Arakawa  and  Sakaki[7].  The  values  at  5  and  10  are  roughly 
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twice  chose  at  2  and  8,  respectively,  which  comes  from  the  double  degeneracy 
of  the  eigenstates. 

V.  Concluding  remarks 

Considering  a  quantum  wire  with  a  rectangular  cross-section  of  a  x  b  and 
very  long  length  L,  which  may  be  called  a  two-direction  DBRTS  when  the 

3 

confining  potential  is  not  very  high,  within  a  very  large  box  of  volume  L  ,  we 
have  calculated  the  local  DOS  and  DOS  in  the  well.  The  latter  shows 
crossovers  from  a  3D  square-root  behavior  to  a  ID  sawtooth-type  behavior,  via 
a  2D  staircase- like  behavior,  when  the  confining  potential  is  very  high  (Uj^  » 
1) .  The  higher  values  of  in  the  case  of  the  transition  from  2D  to  ID 
correspond  to  increased  sharp  peaks  and  finally  reach  the  ideal  sawtooth- type 
behavior  with  singularities  at  the  values  of  E/E^  -  2,  5^  8,  10,  If  we 

3 

consider  a  small  quanuim  box  of  volume  a  x  b  x  c  in  a  large  box  of  volume  L  , 
we  see  a  transition  of  a  finite  DOS  from  ID  to  OD  which  is  expressed  as  a  sum 
of  delta  functions  [7] . 

Although  our  calculations  have  been  performed  so  far  for  rather 
artificial  delta-profiled  systems  only,  we  are  quite  positive  that  this  kind 
of  DOS  transition  will  also  occur  in  real  systems  where,  for  example,  the 
barriers  have  finite  widths.  For  barriers  with  finite  thickness,  the 
effective  mass  of  the  electron  changes  in  passing  from  the  quantiun  well  region 
(GaAs)  to  the  barrier  regions  (AlGaAs)  of  the  structure.  For  this,  BenDaniel 
and  Duke  [8]  suggest  that  current  conservation  is  guaranteed  on  both  sides  by 
use  of  the  boundary  condition 


21  _  21 

m^  dz  m2  dz 


(33) 
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where  and  are  the  effective  masses  of  GaAs  and  AlGaAs,  respectively. 
Also,  Bruno  and  Bahder  [6]  have  considered  this  for  the  one-direction  DBRTS 
case  and  showed  that  the  DOS  at  the  low-energy  subband  edges  is  higher  than 
the  DOS  would  be  at  the  same  energies  in  the  absence  of  barriers  (for  delta- 
profiled  barriers).  So  for  our  two-direction  DBRTS  case,  we  can  estimate  that 
our  result  for  the  DOS  will  also  be  increased  a  bit  upward  at  the  same 
energies  because  of  the  additive  form  of  the  potential  which  we  have  taken. 
Additive  forms  of  potentials  are  used  to  describe  the  motion  of  an  electron  in 
parabolic  quantum  wires  [9]  or  quantum  boxes  [10].  Even  if  we  take  a  quanttim 
wire  with  circular  cross  section,  the  wavefunctions  are  expressed  in  a 
different  way  with  Bessel  functions,  but  the  main  feature  of  our  calculation 
will  not  change  much. 
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Figure  Captions 


1.  Behaviors  of  G  (q)  and  G  (q)  for  U  -  1  and  U  -  10. 

2.  Crossover  of  the  DOS  from  3D  to  2D  in  the  ranee  from  U  -  0  to  U  - 

A  0  0 

2  2 

£ 

20,  i.e.  ,  rrr —  N(E)  as  a  function  of  Here  U  takes  the  values  0,  2, 

zoLm  E.  o 

c  o 


12,  16,  20. 


3.  Crossover  of  the  global  DOS  from  2D  to  ID.  Here  we  take  U-  -  <*>  and  U  - 
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